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The quantum phase transitions in the one-dimensional asymmetric Hubbard model are investi- 
gated with the bosonization approach. The critical conditions for the transition from density wave 
to phase separation, the correlation functions and their exponents are obtained analytically. Our 
results show that the difference between the hopping integrals for up- and down-spin electrons is 
crucial for the happening of the phase separation. When the difference is large enough, the phase 
separation will appear even if the on-site interaction is small. 
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I. INTRODUCTION 

The Hubbard model (HM)i is one of the simplest non- 
trivial models of interacting spin- 1/2 electrons on a lat- 
tice. Its Hamiltonian reads 
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where cj a and Cj tCr , a =|, \ are creation and annihilation 
operators for electrons with spin a at site j respectively, 
n a = c\c a , t is the hoping integral, and U denotes the 
strength of on-site interaction. In one dimension (ID) the 
HM can be solved exactly by the Bethe-ansatz method^ 
The wave function and the energy of N = Nf + N± elec- 
trons on a chain with L sites can be written in terms of 
N pseudo- momentum variables and N± spin rapidities. 
Although the energy spectra have been known for many 
years, the calculation of the correlation functions proved 
to be a delicate problemi 4 -! 5 -! 6 -^ The numerical evaluations 
of the correlation functions 5 and the analytic results 6 - in- 
dicated clearly that the ID HM is a Tomonaga-Luttinger 
liquid (TLL)4£ Assuming that the ID HM is TLL, it 
then becomes possible to calculate the correlation func- 
tions from the knowledge of the energy spectral. Using 
this procedure Schulaii studied the correlation-function 
exponents for different U and band filling n. It is also 
showni^ that the large scale behavior of the spin and 
charge degrees of freedom can be described by two de- 
coupled boson field theories with dynamics governed by 
the TLL Hamiltonian in the small and large U regimes. 

Another nontrivial model is the Falicov-Kimball 
modeli 3 - which consists of localized ions and itinerant 
spinless fermions. The Hamiltonian of the Falicov- 
Kimball model reads 
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where cj are creation operators for spinless fermions, 
and the configuration {wj} denote spatial distribution 



of ions. Clearly, the FKM can be viewed as a modifica- 
tion of the HM in the sense that the one kind of fermions, 
such as down-spin fermions, has infinite mass, and hence 
does not move. Nevertheless, the physics of the FKM 
is completely different. In the neutral case where each 
particle concentration equals 1/2, it was 1 ^ proved that 
the system always orders in an alternating "chessboard" 
phase at a finite transition temperature in all dimensions 
greater than 1. This ordered phase can be interpreted 
as the transition from a high-temperature homogeneous 
(liquid/gas) phase to a low temperature ordered (solid) 
phase. Freericksi 5 - showed that the model (on a hypercu- 
bic lattice) also displayed incommensurate order, segre- 
gation or phase separation. The ID case of the FKM has 
also been extensively studied. Since there is no finite- 
temperature phase transition, the system can have phase 
transition in the ground state. The numerical solutions 16 
produced a conjecture for the case n e + rii < 1 with 
n e = N e /N, rii = Ni/N and the screened Coulomb inter- 
action U is large enough, the system will segregate into 
an empty lattice (with no ions and all the electrons), and 
a full lattice (with all the ions and no electrons). This 
conjecture was later proven to be true by Lemberg erii. 
For any dimensional HFM, Freericks eta6^ gave a theo- 
rem that the strong correlation can leads to PS. 

The relation between the HM and the FKM is straight- 
forward. In order to have a unified framework, the 
asymmetric Hubbard model (AHM) has been introduced 
naturally . 19 ' 20 : 21 ' 22 Its Hamiltonian reads 

L L 

Hahm = ~E E E + C/ E n Jd n id' ( 3 ) 
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where t a is cr-dependent hoping integral. Clearly, if 
f T = t|, the AHM becomes the HM, and if t L = 0, it 
becomes the FKM. The Hamiltonian © has U(l)(g)U(l) 
symmetry for general t a , and the electron number Ax , 
are conserved respectively. In the condensed matter 
physics, the AHM is believed to describe many physical 
phenomena, such as superconductor, valence fluctuating, 
and heavy fermions . 23 ' 24 In the recent development of the 
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optical lattice, it has been pointed out that the AHM can 
be used to describe a mixture of two species of fermionic 
atoms in optical lattices . 25 ' 26 

According to the fact that the ground states of the 
Hamiltonian ([3]) in its two limiting cases: the HM (if = 
t[) and the FKM (if = 0) belong to two different univer- 
sality classes, a quantum phase transition was suggested 
to happen in the phase diagram defined in the U — 
planej 20 ' 22 ' 25 Nevertheless, the quantitative phase dia- 
gram for the case away from half-filling had never been 
obtained until a recent work by Gu, et a! 26 . In their 
work, the quantum entanglement 2 ^ between a local part 
and the rest of the system, and the structure factor of 
charge-density-wave (CDW) for down-spin electrons are 
used to identify the critical point. Here, we are going to 
study the ground-state phase diagram of the AHM away 
from the half filling with the strategy of bosonization 
method. 28 Differ from the numerical approach 26 which 
captures the physics from the finite-size analysis for small 
systems, our work aim to give definitely analytical re- 
sults from the point of view of field theory. The pa- 
per is organized as follows. In section [H] we derive the 
bosonized form of the ID AHM, and clarify the role of 
the some terms in the Hamiltonian. In section IIII1 we 
first diagonalize the effective Hamiltonian in which some 
irrelevant terms are ignored, then obtain the instability 
conditions for the phase separation and compare them 
with the numerical results of a finite sample. We also 
obtain the analytical expressions for the correlation func- 
tions of charge-density- wave, spin-density- wave (SDW), 



singlet-superconductivity (SS), triplet-superconductivity 
(TS) fluctuations, as well as the corresponding expo- 
nents. Finally, a brief summary and acknowledgement 
are given in section [TV] 



II. THE BOSONIZED FORM OF 
ONE-DIMENSIONAL ASYMMETRIC HUBBARD 
MODEL 

The convenient way to analyze the ID AHM is to 
bosonize the Fermi operators and convert them to a quan- 
tum theory of two Bose field a 28 ' 29 . In the framework of 
the standard bosonization method the HM are expressed 
in terms of canonical Bose fields and their dual counter- 
parts as, 
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Here, the Bose field 4> Cl (t>s present the charge and spin 
degrees of freedom, respectively. fc_Ff and kFi are the 
Fermi momentum for up- and down-spin electrons, kp f — 
Tmij/a,kpi — nni/a, with = N^/L and uj, = N^/L 
are the filling densities for up- and down-spin electrons, 
respectively; and a lattice constant. v c jS are the propaga- 
tion velocities of the charge and spin collective modes of 
the decoupled model (Sv = 0), and K c s are the stiffness 



constants. 

It can be seen that this model is the standard HM at 
unpolarized case with both Bose fields decoupling. And 
if model is at half-filling band, kp^ + kp^ = n/a, the 
Umklapp term (the last term in Hb) is important and 
the HM is in the SDW phase with U > or in the CDW 
phase with U < 0. As the model shifts away from the 
half-filling band, the Umklapp interaction can be can- 
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celled. The difference of hopping integral and the filling 
densities of up- and down-spin electrons in the system 
appear as an effect that breaks the spin-charge separa- 
tion as reveals the presence of the third term in the last 
equation. 

Here, we discuss the quantum transition in the case of 
away from the half filling which means that the Umklapp 
interaction can be ignored. In general t^,t± > 0. In the 
parameter region of U > 0, since K s > 1, so the term 
of cos [y^K(j) s j is irrelevant whenever the system is at 
unpolarized case, in the one-loop approximation we can 
cancel it directly and get 



III. RESULTS AND DISCUSSIONS 
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The Hamiltonian (fTTj|) can be diagonalized in terms of 
two new phase fields which contains a mixture of spin and 
charge degrees of freedom. The propagation velocities of 
these collective modes are 
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At these points, the freezing of the lower bosonic mode is accompanied by a divergence in the charge and spin response 
functions. The static charge compressibility k diverges at Sv =Svi, or Sv =Sv2 ■ it behaves as 



As Sv — > 0,u_|_ — > max(v c ,v 8 ) and t>_ — > min (v c , v s ). As Sv increases, w_ decreases until vanishes at the points 
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Beyond these points the susceptibilities becomes negatives. This behavior of the static response functions together 
with the vanishing of the collective modes velocity indicates that the system becomes unstable^ and undergoes a first 
order phase transitional. This instability is known as phase separation and has been shown to occur in the extended 
HM^iM and in the t — J model^. 
In our case, we obtain 



Sv 1 — y v c v s K c K s = at-\ sin(fcpja) + tj sin(fcf^a), 
Sv 2 = 
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It is obvious that Sv\ > Sv%, so the system is in PS phase state as 
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Then we obtain the PS phase state condition Let us now focus our attention on the correlation func- 

, jj\2 tions. Our interest in this work is observing the algebraic 

(^sin(fc n a))(i T sin(A; FT a))< ( — J . (19) 



decay of the instantaneous correlation functions at zero 
temperature and studying how the exponents get mod- 
ified from the standard HM. The operators for CDW, 
SDW, SS, and TS fluctuations in their bosonized form 
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which represent fluctuations of CDW, SDW, SS and TS 
phases, respectively. The correlation functions are de- 
fined as 



Ri(x) = (:Oi(x)Of{0) :). 



(21) 



After some cacluation, we find that the correlation func- 
tions behavior as 
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The exponents a^s determine the divergence of the cor- 
responding phase. The expressions obtained for the a, 
are 
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FIG. 1: (color online) The structure factor of the CDW for 
two different modes q — ir/L and q = N^tt/L at given U /t-\ = 

Q C TJ r in AT AT O + +/+ „„J AT A 4-„~ 



0.05 



0.04 



0.03 



0.02 



0.01 



Bosonalization method 
Exact Diagonalization 




1.5 



FIG. 2: (color online) The boundary line between the den- 
sity wave and phase separation predicted by the bosonization 
method (solid line) and exact diagonalization (dotted line for 
N — 10, JVf = iVj = 2) at give concentration n — 2/5. Here 
t = tJt h u = U/t v 



The ground state is controlled and named by the most 
divergent correlation function, i.e., with the largest on. 

To have a deep understanding on the PS, it is very 
useful to study the structure factor of the CDW. Since the 
dominating configuration of electrons with spin-down is 
quite different in two phase, so we introduce the following 
structure of down-spin electrons, 



'CDW 



(9) = ^E ei9(j "°(Kl"u)-(ni) 2 ) 5 (33) 



where q = 2nir/L, n = 0, 1, • • • , L. We show the struc- 
ture fact as a function of t = tjt-^ at a given U /it = 0.5 
for two different modes q = tt/L and q — Nitt/L in 
Fig. [T] From the figure we find that in the small t 
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FIG. 3: (color online) The boundary line between the den- 
sity wave and phase separation predicted by the bosonization 
method (solid line) and exact diagonalization (dotted line, 
N — 12, iVf = N± = 2) at give concentration n = 1/3. Here 
t = t l /t- t ,u = U/t 1 . 

limit, S(tt/L) dominate, this fact manifest the phase sep- 
aration; while in a relatively larger t region, S(Nin/L) 
dominates, this suggests the density-wave state. So we 
can use the intersection point of the structure factors of 
two different modes to determine the transition point. In 
Fig. [3 and [31 we shown the phase diagram for different 
concentration n — 2/5,1/3 in the small U region. The 
results are very impressive. In both figures, we can see 
that if U < 0.5, the numerical results from the exact 
diagonalization method agree with Eq. (|19[) excellently. 
That is the phase boundary in the small U regime is pro- 
portional to U 2 . However, when U becomes large, say 
U > 1, the bosonization results deviate from the numer- 
ical results apparently. We interpret it due to the fact 
the bosonization method becomes invalid in the large U 
region. On the other hand, the excellent agreement be- 
tween the results obtained from two approaches suggests 
that finite-size correction to the numerical data for a fi- 
nite sample is very small. 



Therefore, the bosonization results are wonderful in 
the small U region and low concentration conditions. 
Since the ID AHM is equivalent to the FKM if t l = (or 
if =0), which has been proven that there exists PS phase 
at infinite-?/ limit when the system shifts away form half- 
filling. From Eq. (|19[). we find that the PS phase always 
appears in the ID FKM whenever the onsite interaction is 
small or large as the system shifts away form half-filling. 
On the other hand, the numerical studies^ suggest there 
might exist a critical U if the density of electrons is close 
to half-filling. This inconsistence may due to the effect 
of Umklapp process around half-filling. 

IV. SUMMARY 

In summary, we have studied the quantum phase tran- 
sitions in the ID AHM with the bosonization approach. 
In the framework of standard bosonization method, we 
first obtained an effective Hamiltonian of ID AHM. Then 
we diagonalized the Hamiltonian and obtained the prop- 
agation velocities of the collective modes for both spin 
and charge degree of freedom. Based on the instability 
condition, we got the final critical conditions of the phase 
transition from DW to PS. We also obtained the analyt- 
ical expressions for the correlation functions of CDW, 
SDW, SS, TS fluctuations, as well as the corresponding 
exponents. 

Our results show that the difference between the hop- 
ping integrals for up- and down-spin electrons is cru- 
cial for the happening of the PS. When the difference 
is large enough, the phase separation will appear even if 
the on-site interaction is small. In the small-?/ and low 
concentration region, critical conditions which scales like 
ti oc U 2 agree with the numerical results excellently. 
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